What Influence Could the Acceptance of Visitors Cause on the Epidemic Dynamics of a Reinfectious Disease?: A Mathematical Model

The globalization in business and tourism becomes crucial more and more for the economical sustainability of local communities. In the presence of an epidemic outbreak, there must be such a decision on the policy by the host community as whether to accept visitors or not, the number of acceptable visitors, or the condition for acceptable visitors. Making use of an SIRI type of mathematical model, we consider the influence of visitors on the spread of a reinfectious disease in a community, especially assuming that a certain proportion of accepted visitors are immune. The reinfectivity of disease here means that the immunity gained by either vaccination or recovery is imperfect. With the mathematical results obtained by our analysis on the model for such an epidemic dynamics of resident and visitor populations, we find that the acceptance of visitors could have a significant influence on the disease’s endemicity in the community, either suppressive or supportive.


Introduction
As the world becomes more of a global village with advances in technology and easier accessibility to different places, it is very crucial to consider side effects like the spread of diseases (Cossar 1994).The history of man is replete with stories of epidemics invading groups of people, sometimes resulting in mortality.In the long 3 Page 2 of 42 run, such diseases can disappear and recur in the future or become less deadly due to people getting immune.Some notable epidemics in history include the "Spanish" flu (1918)(1919), the Black Deaths (1346-1350) which invaded Europe from Asia and recurred for three decades afterwards before getting eliminated (Brauer 2017), the SARS beginning with some infection on an airplane in 2003 (Wang and Wu 2018), and today's pandemic situation with the COVID-19 since mid-December 2019 after the outbreak in China (CDC 2022;ECDC 2022;NIID 2022;WHO 2022).The exponentially increasing number of cases and large-scale spread of the emerging virus about the COVID-19 are being initiated and promoted by the human mobility in global and local scales (Walters et al. 2018;Du Toit 2020;Liu and Saif 2020;Munster et al. 2020;Phan et al. 2020;Ramaswamy et al. 2021;Zhang et al. 2022).
There have been many investigations concerning the effect of a people's displacement due to social and political unrest as well as the natural migration of disease vectors to new areas on the epidemic outbreak, and especially conducted have been many theoretical/mathematical studies taking into account the possibility of individuals becoming infective during transportation and contributing significantly to transport-related infection (see Wilson (2010) and references therein; especially for the SARS virus transmission, see Wang (2014)).Not only the particular transportation with a long travel, but also the human quotidian mobility as a common phase of the human activity can be considered as one of relevant factors that could cause the spread of a transmissible disease such as influenza (WHO 2018;Seno 2020).So is the case of today's pandemic of COVID-19 in local regions of every country (CDC 2022;ECDC 2022;WHO 2022).In the work presented by Parikh et al. (2013), a synthetic population model of the Washington DC metro area was extended to include leisure and business travelers classified as transients.The final size of the epidemic among residents was found to be remarkably higher when transients were included in the simulation of a flu-like disease outbreak.In considering the emerging diseases of wildlife, Tompkins et al. (2015) show that the key drivers of such diseases are agents from domestic sources and human-assisted exposure to infectious agents from wild populations.Talking about swine fever otherwise known as hog cholera, wild boar populations are known to serve as reservoir for the disease, thereby constituting a great challenge for domestic pig farmers, veterinarians and other stakeholders (Mur et al. 2018;Postel et al. 2018).
In this paper, we shall focus on the influence of temporal visitors in a community according to the endemicity of a transmissible disease spreading over the community.Since the globalization in business and tourism becomes crucial more and more for the economical sustainability of local communities, the condition about the acceptance of visitors would be an important part of the community's policy for the public health about a spreading transmissible disease in and out of it.There must be such a decision on the policy by the host community as whether to accept visitors or not, the number of acceptable visitors, and the condition for acceptable visitors.Actually the importance of such a policy on the tourism regulation has been recognized more and more in the post-COVID-19 period (Rastegar et al. 2021;Volgger et al. 2021;Yan et al. 2021;Jones 2022;Okafor and Yan 2022).Making use of a mathematical model modifying the basic SIR model incorporating the regulated acceptance of an amount of temporal visitors in the community, which may be called an SIRI type of mathematical model, we shall try to consider the influence of visitors on the spread of a reinfectious disease in a community, especially assuming that a certain proportion of visitors are immune at the entry into the community.
The reinfectivity of disease in this paper means that the immunity gained by either vaccination or recovery is imperfect.For a spreading transmissible disease accompanied with a reinfectivity, the acceptance of visitors must influence the endemicity of such a disease in the community.Then the community's policy must take account of the reinfection risk for both of residents and visitors.Actually there are transmissible diseases with a reinfectivity, including influenza (Davies et al. 1984;Hay et al. 2001;Earn et al. 2002;Price et al. 2022;Wang et al. 2022), pertussis (Hethcote 1999;van Boven et al. 2000), Lyme disease (Nadelman et al. 2012), hand, foot and mouth disease (Zhang et al. 2019), malaria (Arias et al. 2022;Rehman et al. 2022), tuberculosis (Vynnycky and Fine 1997;Horsburgh et al. 2022;Qiu et al. 2022), Ebola virus disease (MacIntyre and Chughtai 2016; Agusto 2017), chronic lung diseases (Yum et al. 2014), invasive pneumococcal disease (Lipsitch 1997), meningococcal disease (Gupta and Maiden 2001), and COVID-19 (Crawford 2022;Kumar et al. 2020;Le Page 2022;Mensah et al. 2022;Nguyen et al. 2022;Ren et al. 2022;Saad-Roy et al. 2022;Salzer et al. 2022;Shaheen et al. 2022), although the reinfectivity has been still requiring scientific researches to understand its kinetics and other nature.
According to Chowell et al. (2016), it is crucial to formulate reliable models that embody the basic transmission characteristics of specific pathogens and social scenarios.They further stated that improved models are required to capture the variation in early growth dynamics of real epidemics in order to gain better understanding of the dynamics as they reviewed trends in modeling and classifying early epidemic progression.Recently the mathematical epidemic dynamics models are being used to estimate or evaluate some epidemiological parameters and to predict the temporal variation in the morbidity about a spreading disease, making use of epidemiological data (Siettos and Russo 2013).So are particularly those on the COVID-19 spread (for example, Athayde and Alencar 2022;Kobayashi et al. 2020;Lin et al. 2022;Musa et al. 2022 ), however, this is not the case in our paper.
We are going to try to shed a light on the theoretical side about the influence of temporal visitors on the epidemic dynamics with a transmissible disease with a reinfectivity, since the acceptance of visitors under such an epidemic dynamics with the risk of reinfection could be a crucial factor to influence the endemicity in the community which must take a policy about the acceptance (Dansu and Seno 2019;Crawford 2022;Salzer et al. 2022).We shall consider one of the simplest epidemic dynamics models for a transmissible disease with reinfectivity, which is an SIRI type of mathematical model, taking into account the influence of temporal visitors accepted by a community.With the mathematical results obtained by our analysis on the model, it will be implied that the acceptance of visitors could have a significant influence on the disease's endemicity in the community, either suppressive or supportive, depending on the risk of reinfection and the nature of accepted visitors.

Assumptions
We consider the spread of a transmissible disease during a short-term period, that is, a season after the community starts to accept visitors from the outside, satisfying the following assumptions on the epidemic dynamics:

[H1]
The demographic change in the resident population is negligible in the season. [H2] The fatality of disease on the resident and visitor populations is negligible in the season. [H3] The community starts the acceptance of a number of temporal visitors from the outside in the season after a transmissible disease has already invaded in it. [H4] The entry flow of visitors is constant, that is, the net entry rate is constant independently of time. [H5] The exit of visitors from the community follows a constant per capita exit rate.

[H6]
No infected visitor is accepted by the community (i.e., the perfect quarantine), so that every accepted visitor is susceptible or immune to the disease at the entry into the community.

[H7]
A given proportion of visitors is immune at the entry into the community.

[H8]
Only the susceptible residents can get the vaccination to become immune, and it is not available for any visitor staying in the community.

[H9]
Immune visitor has a possibility to get reinfected (i.e., the imperfect or partial immunity) during its stay in the community, the same as the immune resident does.
[H10] Infected visitor has the same exit rate as the susceptible visitor, that is, we neglect any influence of the infection on the visitor's stay in the community.
Assumption H1 indicates a time-independent constant size of resident population during the season in which the epidemic dynamics is going on.We then ignore the death due to the transmissible disease under consideration in the epidemic dynamics too, as indicated by the assumption H2.Assumption H3 indicates that the community accepts the visitors, even undergoing the spread of a transmissible disease, since the fatality of the disease is negligible with the assumption H2.No disease invasion with the visitors is assumed, as indicated by the assumption H6.From the assumptions H4 and H6, the community carries out the perfect regulation for the visitors at the entry concerning the entry number and the quarantine.Assumption H5 mathematically means that the exit of a visitor from the community follows the homogeneous Poisson process.In a model with ordinary differential equations, it can be introduced with a constant exit rate per visitor.Assumption H7 is to reflect the situation of public health out of the community, applying the mean-field approximation for the proportion of immune visitors at the entry.Since we assume that the community undergoes the disease spread, the assumption H8 gives the existence of a vaccination program for the residents, while it is not applied to the visitors.However, since the disease is reinfectious as assumed by the assumption H9, the immunity obtained by the vaccination or the recovery from the disease works only to reduce the risk of reinfection.Hence the state transition in terms of the disease follows the susceptible-infective-recovered/immunized-infective (SIRI) structure in our modeling, as used for example in Gomes et al. (2004Gomes et al. ( , 2005)), Gökaydin et al. (2007), Stollenwerk et al. (2007), Martins et al. (2009), Pinto et al. (2010), Song et al. (2011), Georgescu andZhang (2013), Guo et al. (2014), Pagliara et al. (2018), Buonomo (2020), Ghosh et al. (2020), Wang (2021), Srivastava et al. (2022).
Remark that the assumed reinfection is not caused by the waning or loss of immunity, which must take a certain period after getting it by the infection or vaccination.As already mentioned in the introduction section, we assume instead the imperfectness of immunity obtained by the infection or vaccination.Hence we do not introduce any specific period or time scale to get reinfected after getting the immunity.Since the infection or vaccination must generate an immunity against the disease, the assumption H9 indicates that the immunity is imperfect or partial against the infection, for example, due to the multiplicity of pathogen types (e.g., mutated variants) (Gökaydin et al. 2007;Wang et al. 2022).Because the cross-immunity is wellknown for such similar pathogens by the antigen for a type of pathogen, the reinfection may be suppressed or fail to induce an effective symptom to reproduce and discharge the pathogen to the environment.
For a simplification, the assumption H10 indicates that the exit of visitor is independent of whether the visitor is infected or not during the stay in the community.This assumption would be appropriate when the expected duration of the visitor is sufficiently shorter than the latent period, whereas it may be less appropriate when it is long.As assumed by the assumption H2, we consider a transmissible disease with little serious symptom, so that the assumption H10 would be applicable for visitors infected by such a disease.

Generic model
With these assumptions given in the previous section, we shall consider the following model of ordinary differential equations (Fig. 1): where S v , I v , and R v are the subpopulation sizes of susceptible, infective, and immune visitors respectively.Similarly S r , I r , and R r are the corresponding subpopu- lation sizes about the residents.The population sizes of residents and visitors staying in the community are denoted by N = S r + I r + R r and m = S v + I v + R v respec- tively.The resident population size N is constant independently of time t, as seen from d(S r + I r + R r )∕dt = 0 for any t by the system (2).From the assumption H3 in Sect.2, the community starts the acceptance of visitors from the outside in the considered season after a transmissible disease has already invaded in it.The visitor population size m could be reasonably assumed to be less than the population size of residents N: m < N , whereas we shall not specifically assume so but consider the mathematically general case of m in the subsequent sections without any constraint except for m ≥ 0 .On the other hand, as given in Sect.3.2, we will take an assump- tion on the visitor population size m accompanying with a confinement for the net entry rate of visitors Λ.
All parameters are positive.Parameter is the proportion of immune visitors at the entry ( 0 ≤ ≤ 1 ).Proportion 1 − of visitors is susceptible at the entry.Param- eter q is the per capita exit rate of visitor.Thus the expected duration of a visitor's stay in the community is given by 1/q.Parameter is the reinfection coefficient for immune resident and visitor, while is the infection coefficient for susceptible ones.Then the infection forces for the susceptible individual and the immune individual are respectively given by (I r + I v )∕(N + m) and by (I r + I v )∕(N + m) for both resident and visitor.Remark that, in the setup for our modeling, the visitors do not form any specific subcommunity distinct from the resident population.From the assumption H3, they are temporal visitors for tourism, business etc.For this setup, we could assume that most of visitors are independent of the others.Thus, for a mathematical simplification, the influence of their movement on the epidemic (1) Dynamics for the visitor population: (2) Dynamics for the resident population: dynamics is introduced in the epidemic dynamics by the mean-field approximation.Further, although the visitors' contribution to the infection forces would be different from the residents' one because of the difference in the mobility/behevioral pattern, the infection forces have their same contributions in our modeling here.This modeling may be regarded as an oversimplification, though we think that our model would still worth being considered to get a cue for the discussion about the influence of visitors on the epidemic dynamics within a community.
From the assumption H9 in Sect.2, our modeling assumes that the immunity is imperfect or partial against the infection.Because of the cross-immunity by the obtained antigen, we reasonably assume that 0 ≤ ≤ 1 in our model, so that the reinfection coefficient is not beyond the coefficient for susceptible .That is, the infection after the vaccination or recovery from the disease generally has a smaller likelihood than that for the susceptible.For the extremal case of = 1 , the vaccination or recovery does not work at all to reduce the risk of reinfection.For = 0 , the recovery and vaccination give the perfect immunity so that there is no likelihood of reinfection.Thus the parameter means an index for the likelihood of reinfection after the recovery or vaccination, so that it can be regarded as an index for the risk of reinfection.Remark here again that the reinfection in our modeling is assumed to be not due to the waning of immunity (like for the SIRS models) but due to the imperfect immunity, and hence also the vaccinated individual has a risk to get infected, as introduced by the assumption H9.
Parameter is the recovery rate of an infective individual, and the recovered individual gets immunity, which is however imperfect.Only the susceptible residents can get the vaccination, with rate , and it is not available for any visitor staying in the community.Since the vaccination is imperfect from the assumption H9, it works to reduce the risk of infection but is unable to protect the vaccinated individual from the infection.

Assumption for the visitor population size in the community
According to the dynamics for the visitor population (1), we have is the visitor population size at time t, Λ the net entry rate of visitors, and q the per capita exit rate of visitor.Now let us consider the stationary situation with respect to the temporal change of visitor population size.This means that the number of visitors is assumed to be stationary, which may be regarded as a consequence of the regulation of their entry by the community, following the assumptions H4 and H6 given in Sect. 2. Therefore we mathematically assume the situation to satisfy that dm∕dt = 0 .Hence we put Λ = qm , and hereafter treat the visitor population size m as a positive constant.

The Initial Condition
Following the assumption of the stationary visitor population size with Λ = qm as given in Sect.3.2, we have the following dynamics for the visitor population at the disease-free state: where R v means the subpopulation size of immune visitors in the disease-free com- munity.It can be easily found that this dynamics results in an eventual approach to the equilibrium state such that For this reason, let us assume the following initial condition for the epidemic dynamics with the model given by the system of ( 1) and ( 2): where This initial condition defines the situation when the community starts the acceptance of visitors from the outside, even under the existence of disease in it.The setup of this initial condition as our modeling follows the assumption H3 in Sect. 2.

Basic Reproduction Number
In the biological/epidemiological context, the basic reproduction number is defined as the expected number of new cases of infection caused by an infective individual within a population consisting of susceptible contacts only (for a useful review about the definition, the translation, and the practical application, see Delamater et al. (2019)).Such a situation can never occur for the epidemic dynamics in reality or even by any mathematical model.This is because the secondary infection itself changes the situation of population where the disease spreads.Increase of infectives must reduce the likelihood of the contact between an infective and an susceptible since the likelihood to contact the other infectives becomes less negligible.For this reason, the basic reproduction number must be mathematically defined to match the above conceptual definition in the biological/epidemiological context.Briefly saying, it is mathematically defined as the supremum of the expected number of new cases of infection caused by an infected individual at the stage of disease invasion in the community (Seno 2022).Following this definition in a biological/epidemiological sense, a mathematical theory can be used to derive the basic reproduction number about the mathematical model of epidemic dynamics, for example, as the spectral radius of a specific matrix, called the "next generation matrix", for a system of ordinary differential equations governing an epidemic dynamics (see Diekmann et al. (2013) for a complete reference, or see Lewis et al. (2019) or van den Driessche (2017) for the further review).The choice of the derivation way of R 0 may depend on the nature of the model.As well known on the nature of the basic reproduction number, the invasion of a transmissible disease in the community is successful if R 0 > 1 , while it fails if R 0 < 1 .In this context, the invasion success means the increase of infectives in the early stage of the disease spread after a sufficiently small number of infectives appear in the community, and the invasion failure implies the decrease.
For our model given by the system of ( 1) and ( 2), we can derive the following formula of the basic reproduction number R 0 (Appendix A): where ∶= m∕N , and for a convenience in the following arguments, we define R 00 ∶= ∕ , which is the basic reproduction number for the community when no visitor comes in (i.e., m = 0 ).Note that this basic reproduction number is funda- mental for the epidemic dynamics in the community after it starts the acceptance of visitors.
From the formula (4), we can immediately find that the basic reproduction number R 0 is less than 1 independently of the nature of accepted visitors if R 00 < 1 .Hence, if the disease fails its invasion in the community with R 00 < 1 before starting the acceptance of visitors, the number of infectives in the community cannot turn to increase in the early period after the acceptance of visitors starts.As we will see in the later sections of the analysis on our model, this is valid only in the early period after the acceptance of visitors starts.
As for the dependence of R 0 on the nature of accepted visitors, we note that R 0 is monotonically decreasing in terms of when the visitors contain some immune ones (i.e.,  > 0 ).Moreover, R 0 becomes smaller as the proportion of immune visitors at the entry gets larger.If any visitor is susceptible, that is, when = 0 , there is no contribution of the visitors to the basic reproduction number R 0 , that is, R 0 = R 00 .
As an extremal case, if the immunity gained by the vaccination or recovery from the disease does not work at all to reduce the risk of reinfection, that is, if = 1 , the basic reproduction number R 0 is independent of the acceptance of visitors.This is easily understandable, since the reinfection is regarded as the same as the infection for the susceptible so that the immune individual is regarded as equivalent to the susceptible according to the epidemic dynamics when = 1 .Such an extreme case may be regarded as corresponding to an SIS type of the epidemic dynamics, where the state transition in terms of the disease follows the susceptible-infective-susceptible structure. (4) the supremum of the expected new cases per unit time for the resident.
the supremum of the expected new cases per unit time for the visitor, given by the sum of secondary infections for susceptible and immune visitors.
As the other extremal case, if the immunity is perfectly effective to make the immune individual unable to be reinfected, that is, if = 0 , the entry of immune visi- tors works to reduce the value of R 0 for the community.This extremal case may be regarded as corresponding to an SIR type of the epidemic dynamics, where the state transition in terms of the disease follows the susceptible-infective-removed structure.
From these arguments with the basic reproduction number R 0 given by ( 4), we can get the following result on the influence of the acceptance of visitors in the early period after the community starts the acceptance of visitors (Fig. 2): Theorem 4.1 The acceptance of visitors influences the basic reproduction number R 0 given by (4) for the epidemic dynamics with the system of (1) and (2) as follows: (i) The acceptance of visitors makes R 0 smaller than R 00 if and only if the visi- tors contain some immune, and its decline becomes greater as the number of accepted visitors gets larger; ∞ , the acceptance of visitors so many as makes R 0 < 1.
The result (i) in Theorem 4.1 means that the acceptance of visitors does not help the spread of a transmissible disease as long as R 00 ≤ 1 , and instead it could work to sup- press the spread if the community accepts a sufficiently large number of visitors with a sufficiently large proportion of immune, as indicated by the results (ii) and (iii).
Fig. 2 The dependence of the basic reproduction number R 0 given by (4) on parameters (1 − ) and ∶= m∕N .Numerically drawn with R 00 = 1.4 Note that these arguments and result are about the effect of the acceptance of visitors on the temporal change of the number of infectives only in the early period after the community starts the acceptance of visitors.They cannot be necessarily applied for its later temporal change.It may be possible that the number of infectives turns to increase in a later period, independently of what final state the epidemic dynamics approaches, as we will actually see in the later sections of the analysis on our model.

Non-Dimensional Transformation of the System
Since the population sizes of visitors and residents are assumed constant independently of time, the above six dimensional system of ( 1) and ( 2) can be mathematically reduced to the following closed four dimensional one, making use of Now we apply the following transformation of variables and parameters for this four dimensional system: and then we can derive the following non-dimensinalized system: where R 00 ∶= ∕ as before.Remark that the symbol R 00 is formally introduced now as a dimensionless parameter for the non-dimensionalized system given by ( 5), while its meaning is given in Sect. 4 as the basic reproduction number for the community when no visitor comes in.The initial condition (3) now becomes with x r0 ∶= S r0 ∕N = 1 − I r0 ∕N = 1 − y r0 , where we set R r0 = 0 for a simplicity that there is no immune resident at the initial, which will not affect any mathematical result obtained in our analysis on the model.In the following sections, we shall analyze the non-dimensionalized system (5) to investigate the nature of the epidemic dynamics by our model with the system of ( 1) and (2).

Dynamics Without Reinfection
In this section, we consider the system without reinfection, that is, with = 0 , while we will consider our model only with  > 0 in the subsequent sections.For the system (5) with = 0 , we can easily find that x r → 0 and y r → 0 as → ∞ .In other words, since the epidemic dynamics for the resident population is governed by an SIR model with the continuous vaccination measure for the susceptibles, the disease eventually disappears in the resident population, and the residents come to make no contribution to the epidemic dynamics.This means that, for our model of ( 1) and ( 2) without reinfection, no endemic state can be established as long as the community does not accept any visitor from the outside.Thus we consider the case of  > 0 hereafter in this section, when the community accepts the visitors.By the local stability analysis with the eigenvalues of Jacobi matrix at the equilibrium, we can easily find that the endemic equilibrium E +0 (x * v , ỹ * v , 0, 0) with is locally asymptotically stable when it exists.Then we can get the following theorem on the epidemic dynamics given by ( 5) with = 0 (Appendix B): Theorem 6.1 For the system (5) with = 0 , if and only if the condition is satisfied, the endemic equilibrium E +0 (x * v , ỹ * v , 0, 0) with (6) uniquely exists, and it is globally asymptotically stable.Otherwise, the disease-eliminated equilibrium E 00 (1 − , 0, 0, 0) is globally asymptotically stable.
See Fig. 3(a, b) for numerical examples.
From the basic reproduction number R 0 given by (4) with = 0 , we can easily find that R 0 > 1 if the condition (7) is satisfied, while the inverse does not neces- sarily hold.Hence we can get the following result: Corollary 6.1.1Even when the disease successfully invades in the community with R 0 > 1 , the disease without reinfection eventually gets eliminated unless the condi- tion (7) is satisfied.
When the condition ( 7) is unsatisfied with R 0 > 1 , the number of infectives increases at the initial stage of disease spread in the community, and then it eventually turns to decrease toward zero, as numerically exemplified in Fig. 3(a).
For = 1 , the condition (7) cannot be satisfied: Corollary 6.1.2If the community accepts only immune visitors, the epidemic dynamics without reinfection necessarily approaches the disease-eliminated equilibrium.
Therefore, when the reinfection is impossible/negligible, the acceptance of visitors with a high immune proportion at the entry does not cause the endemicity of disease.We must remark that the endemic equilibrium E +0 is sustained only by the visitor subpopulation, while no resident contributes to the epidemic dynamics at the equilibrium because all residents have eventually become immune by the past infection or vaccination [Fig. 3(a,b)].From Theorem 6.1, for the disease with a sufficiently high infectivity, the acceptance of many visitors with a sufficiently small immune proportion at the entry can induce such an apparent endemicity in the community.Therefore, as seen in Fig. 3(a, b), when the community successfully controls and reduces the number of visitors to make sufficiently small, or if the community suspends accepting the visitors, the endemic state can be disrupted, and then the disease gets eliminated in the community.However, as we will see in the subsequent sections, if the disease is accompanied by a reinfectivity, this could not be the case [Fig.3(c, d)].

Dynamics with no Visitor
Next we consider the model of ( 1) and ( 2) with reinfection, that is, with  > 0 , when the community does not accept any visitor from the outside.Thus we analyze the following system derived from (5) with x v = y v ≡ 0 and = 0: As already mentioned in Sect.4, the basic reproduction number for this epidemic dynamics is defined as R 0 = R 00 = ∕ .It is easy to find that there are two feasible equilibria for this system of (x r , y r ) : the disease-eliminated equilibtium E 0 (0, 0) and the endemic equilibrium E + 0, 1 − 1∕( R 00 ) .The endemic equilibrium E + exists when and only when R 00 > 1.
Making use of the local stability analysis for the equilibrium of the system (8), we can easily find that the equilibrium E 0 is locally asymptotically stable if R 00 < 1 , and unstable if R 00 > 1 .The endemic equilibrium E + is locally asymptotically (8) Fig. 4 Application of the isocline method for the system with no visitor (8) when a R 00 < 1 ; b R 00 = 1 ; (c) R 00 > 1 stable whenever it exists.Taking account of the result on the local stability of equilibria, the isocline method for the two dimensional system (8) can further give the following result (see Fig. 4): Theorem 7.1 For the system (8) with no visitor, (i) if and only if R 00 ≤ 1 , the disease-eliminated equilibrium E 0 is globally asymptotically stable; (ii) if and only if R 00 > 1 , the endemic equilibrium E + exists and is globally asymptotically stable, while E 0 is unstable.
This result was shown for a mathematically equivalent SIRI model in Gomes et al. (2005).
From this result, we find that, even with the basic reproduction number R 00 > 1 , the community approaches the disease-eliminated equilibrium E 0 if R 00 ≤ 1∕ with 0 <  < 1 .In such a case, the number of infectives increases at the initial stage of disease spread in the community, and then it eventually turns to decrease toward zero.If and only if the basic reproduction number is sufficiently large as R 00 > 1∕ , the disease becomes endemic in the community.
As shown in the previous section, no endemicity arises in the community with a non-reinfectious disease when no visitor is accepted.Now the result obtained in this section indicates that the endemicity of a disease can arise in the community even with no visitor if the disease has both sufficiently high infectivity ( R 00 > 1 ) and sufficiently high reinfectivity (  > 1∕R 00 ).

Disease-Eliminated Equilibrium
For the system (5) with visitors, we can get the following result on the local stability of the disease-eliminated equilibrium E 00 (1 − , 0, 0, 0) (Appendix C): Theorem 8.1 The disease-eliminated equilibrium E 00 is unstable if while it is locally asymptotically stable if the inverse inequality of (9) is satisfied.
It can be easily seen that the condition (9) becomes equivalent to (7) as → 0 .This shows a mathematical consistency of Theorem 8.1 to Theorem 6.1.
Moreover we can prove the following result about the relation of the basic reproduction number R 0 defined by (4) to the condition (9) (Appendix D): This result is consistent with the biological/epidemiological meaning of the basic reproduction number with respect to the invasion success of a disease in a population, which was described in Sect. 4 and references therein.Under the condition that R 0 < 1 at the initial stage of a disease invasion in a population, the number of infectives is expected to decrease toward the disease elimination.In a sense of epidemic dynamics, such a decline of the infective subpopulation toward the elimination must follow the locally asymptotic stability of the disease-free equilibrium (as referred in most literatures), which corresponds here to the disease-eliminated equilibrium E 00 .The result of Corollary 8.1.1 shows this consistency.

Endemic Equilibrium
From the equations of ( 5), if the endemic equilibrium E ++ (x * v , y * v , x * r , y * r ) with y * v > 0 and y * r > 0 exists for  < 1 , then it must satisfy that and x * r = 0 .In contrast, at the endemic equilibrium E ++ for = 1 , we have and x * v = 0 , x * r = 0 , instead of (10).We can obtain the following result on the exist- ence of E ++ (Appendix E): Theorem 8.2 The endemic equilibrium E ++ uniquely exists if and only if the condi- tion (9) is satisfied.
Hence, when the disease-eliminated equilibrium E 00 is locally asymptotically stable, the endemic equilibrium E ++ does not exist, and when E 00 is unstable, E ++ uniquely exists.
For the local stability of the endemic equilibrium E ++ for (5), we can get the fol- lowing result: (10) Theorem 8.3 When the endemic equilibrium E ++ exists, it is locally asymptotically stable.
This theorem can be proved by the eigenvalue analysis on the Jacobi matrix for (5) at the endemic equilibrium E ++ , applying the Routh-Hurwitz criterion (Appendix F).Although we could not get any analytical result on the global stability of the endemic equilibrium E ++ , our numerical calculations imply that it is globally asymptotically stable when it exists.We then have the mathematical consistency of Theorems 8.1, 8.2, and 8.3 to Theorem 6.1 as for the case with no reinfection, = 0.

Shift in Endemicity
As the important preliminary for our analysis on the model, we can easily find the following features of G( , ) defined in (9): monotonically decreasing in terms of  if  <  s .• G(, ) < 1 for any positive and  <  s if and only if Then from these mathematical features of G( , ) , and Theorems 8.1, 8.2, and 8.3, we can get the following result on the disease endemicity in the community accepting the visitors (see Figs. 5 and 6): Theorem 9.1 Independently of whether the community accepts the visitors or not, it approaches an endemic equilibrium if R 00 ≥ max 1, G ∞ ( ) , while it approaches the disease-eliminated equilibrium if R 00 ≤ min G ∞ ( ), 1 , where For the region Ω + , the acceptance of visitors may change the endemic situation of the community for the disease-eliminated equilibrium as described in Theorem 9.2, while, for the region Ω − , it may drive the situation of the community approaching the disease-eliminated equilibrium toward the endemic equilibrium as described in Theorem 9.3.For the region out of Ω − and Ω + , the endemicity is independent of whether the community accepts visitors or not Only when the endemicity could significantly depend on the acceptance of visitors.
The inequality (12) gives a necessary condition for which the acceptance of visitors could change the epidemic situation in the community from the endemic equilibrium to the disease-eliminated equilibrium or vice versa.The corresponding parameter regions are shown as Ω − and Ω + in Fig. 6.Numerical examples of such a change of endemicity by the acceptance of visitors are given in Fig. 7(a, c).
Further from the monotonicity of G( , ) in terms of and as described in the above, we find the following result on the condition with respect to the influence of the acceptance of visitors on the endemicity in the community:
Moreover, from the features of G( , ) described in the above, we find that, if R 00 ≥ G ∞ ( ) ≥ 1 , the disease is endemic independently of how many visitors the community accepts even under the condition ( 13), as seen in Fig. 5. Thus, in comparison to the result for the community with no visitor (i.e., = 0 ) given by Theorem 7.1, we can get the following result (see Fig. 6): Theorem 9.2 Suppose that the disease was endemic under the condition that 1 < R 00 < 1 + c before the community accepts visitors.If the community accepts visitors with an immune proportion the disease remains endemic independently of how many visitors are accepted.If the community accepts visitors with an immune proportion  >  ∞ , then the acceptance of visitors so many as makes the community approach the disease-eliminated equilibrium.Even if  >  ∞ , the acceptance of visitors with ≤ c does not sufficiently shift the endemicity, and the disease remains endemic.
The critical value ∞ satisfies the equation R 00 = G ∞ ( ∞ ) .The latter case of  >  ∞ in Theorem 9.2 corresponds to the parameter region Ω + in Fig. 6. Figure 8(b) shows a numerical example of the ( , )-dependence in such a case when R 00 > 1.Since c defined by ( 15) is monotonically decreasing in terms of when R 00 > 1 , we note that where  c1 > 0 for 1 < R 00 < 1 + c .Then we get the following corollary: Corollary 9.2.1 When 1 < R 00 < 1 + c , if the community accepts visitors few enough to have  <  c1 , the disease remains endemic independently of how much proportion of visitors is immune at the entry.
This result is numerically pointed out in Fig. 8(b, c).Moreover we note that, when R 00 < 1 + c so that ∞ defined by ( 14) is nega- tive, the first case in Theorem 9.2 does not occur.Then the ( , )-dependence becomes as shown by Fig. 8(c), where there is a finite value of beyond which the community approaches the disease-eliminated equilibrium, independently of the immune proportion in the visitors at the entry: Corollary 9.2.2When 1 < R 00 < (1 + c) , if the community accepts visitors so many as  >  c0 , the disease tends to eliminated independently of how much proportion of visitors is immune at the entry.When (1 + c) ≤ R 00 < 1 + c , only the acceptance of visitors with  >  ∞ and  >  c can change the endemicity and lead the community to the disease-eliminated equilibrium.
Hence the value c0 gives a sufficient number of accepted visitors which is effective to suppress the disease spread in the community when 1 < R 00 < (1 + c) .See the numerical examples in Fig. 7(b, c).
In contrast, when the risk of reinfection is so weak as  < 1∕(1 + c) , the acceptance of visitors may cause the opposite influence on the epidemic dynamics, as numerically indicated by Fig. 8(a): Theorem 9.3 Suppose that the disease was getting eliminated under the condition that (1 + c) < R 00 < 1 before the community accepts visitors.If the community accepts visitors with an immune proportion ≥ ∞ , the disease keeps getting elimi- nated independently of how many visitors are accepted.If the community accepts visitors with an immune proportion  <  ∞ , then the acceptance of visitors so many as  >  c induces the endemicity, and the disease becomes endemic.Even if  <  ∞ , the acceptance of visitors so few as ≤ c does not induce the endemicity, and the disease keeps getting eliminated. .
3 Page 22 of 42 The situation considered in this theorem corresponds to the parameter region Ω − in Fig. 6, and is numerically exemplified by Fig. 7(a).Theorem 9.3 indicates that, if the proportion of immune visitors is so low as  <  ∞ , there exists the upper threshold c for the number of accepted visitors to suppress the revival of the disease spread after starting the acceptance of visitors in the community where the disease was getting eliminated.

Acceptance of Only Immune Visitors
When the community accepts only immune visitors, that is, when  > 0 with = 1 , G( , 1) is necessarily greater than 1 and monotonically increasing in terms of .Then, from Theorem 9.2, we can find that only the acceptance of visitors so many as  >  c1 can induce the disease-eliminated equilibrium in the community where the disease was endemic before starting the acceptance of visitors.As defined by ( 16), the critical value c1 depends on the risk of reinfection, and then, from Theorems 9.1 and 9.2, we can obtain the following result (see Figs. 8 and 9): Corollary 9.2.3 Suppose that the disease was endemic under the condition that R 00 > 1 before the community accepts only immune visitors.If R 00 ≥ 1 + c , the Fig. 9 Parameter region and boundary indicated by the condition (9) with the results given by Theorems 8.1, 8.2, and 8.3 when the community accepts only immune visitors with = 1 .Numerically drawn with R 00 = 4.0 and c = 1.0 .Refer to Sect.9.2 endemicity remains independently of how many visitors the community accepts.If 1 < R 00 < 1 + c , the acceptance of visitors so many as  >  c1 is effective to make the disease eliminated.
Therefore the community under an endemic situation could have a preferable influence to suppress the endemicity by accepting only immune visitors only when the reinfectivity is sufficiently low as indicated by Fig. 9.
If the community was approaching the disease-eliminated equilibrium with the risk of reinfection so low as R 00 ≤ 1 , the community keeps approaching the dis- ease-eliminated equilibrium even after starting the acceptance of only immune visitors, independently of how many visitors the community accepts.

Acceptance of Only Susceptible Visitors
Now let us consider the case where all visitors accepted by the community are susceptible, that is, when  > 0 with = 0 .Then, from Theorem 9.3, we find three different cases according to the influence of the acceptance of visitors as shown in Fig. 10, taking account of the features of G( , 0) given in Sect. 9.1.We can get the following result: Corollary 9.3.1 When the disease was getting eliminated under the condition that R 00 ≤ 1 , the acceptance of only susceptible visitors induces In contrast, when the disease was endemic under the condition that R 00 > 1 , the acceptance of only susceptible visitors induces no endemicity if R 00 ≤ (1 + c); no endemicity if R 00 > (1 + c) and  ≤  c0 ; the endemicity if R 00 > (1 + c) and  >  c0 .The critical value c0 is defined by ( 17).Therefore the acceptance of only susceptible visitors could have the counter effect according to the endemicity, depending on the infectivity of disease.Only for a moderately high infectious disease such that 1∕ < R 00 < 1 + c , the accept- ance of only susceptible visitors so many as  >  c0 can lead the community to the disease-eliminated equilibrium.For the disease with a low reinfectivity such that 1 + c < R 00 ≤ 1∕ , the acceptance of only susceptible visitors so many as  >  c0 can lead the community to the endemic equilibrium.
Further, as indicated by Fig. 10, we find that there is a sufficient value of which determines the epidemic situation in the community: Corollary 9.3.2If R 00 > 1 + c , the acceptance of only susceptible visitors so many as necessarily makes the disease endemic.In contrast, if R 00 < 1 + c , the acceptance of only susceptible visitors so many as necessarily makes the disease eliminated.
The former case means an unpreferable influence of the sufficiently large number of visitors for the community with the spread of a highly infectious disease, while the latter does a preferable influence for the community with the spread of a moderately infectious disease.

Change in Endemic Size
Figure 11 shows the numerically drawn -dependence of endemic sizes y * r and y * v at the endemic equilibrium determined by (10).As the figure implies, the endemic size necessarily has a monotonic dependence on the number of accepted visitors, represented now by , about which we can get the following analytical result (Appendix G): Theorem 9.4 The endemic sizes y * r , y * v , and the total endemic size no change in the endemicity if R 00 ≥ (1 + c); no change in the endemicity if R 00 < (1 + c) and  ≤  c0 ; the elimination of disease if R 00 < (1 + c) and  >  c0 .
are monotonically increasing in terms of if and only if R 00 ≤ 1 or  It can be easily found that  c < 1 when R 00 > 1 .See the numerically drawn ( , )-dependence of the endemic size y * r in Fig. 12.We remark that, as shown in Sect.9.2, the endemic equilibrium for = 1 exists only when R 00 > 1 .Then we can find the following result too (Appendix H): Corollary 9.4.1 When the community accepts only immune visitors (i.e., = 1 ), the endemic size is monotonically decreasing in terms of .This result may be regarded as included in Theorem 9.4 because the condition given in Theorem 9.4 can never hold when = 1 .We can see the numerical exam- ples in Fig. 12(b, c).
When c defined by ( 19) is non-positive with R 00 > 1 , any cannot be smaller than c , so that the endemic size is necessarily monotonically decreasing in terms of : Corollary 9.4.2If R 00 > max 1, 1∕ − c , the endemic size is necessarily mono- tonically decreasing in terms of .
The numerical example of Fig. 12(c) illustrates the case.
For the critical case of  =  c > 0 with R 00 > 1 , we can derive the explicit val- ues at the endemic equilibrium E 00 from (10) (Appendix G): Hence the endemic sizes are independent of the number of accepted visitors in this case: Corollary 9.4.3For  =  c > 0 with R 00 > 1 , the endemic sizes y * r , y * v , and the total endemic size z * are determined independently of .
Additionally we can find the following relations among the specific values s , ∞ , and c for the immune proportion of accepted visitors at the entry: Consequently as indicated by Fig. 13, the larger number of accepted visitors makes the endemic size bigger only when the immune proportion of accepted visitors at the entry is sufficiently small under the epidemic situation with a sufficiently low risk of infection.

Concluding Remarks
The results of our model imply that the acceptance of temporal visitors from the outside may induce a significant change of the epidemic state in the community.Contrary to an intuitive expectation, the acceptance of visitors does not necessarily make the epidemic situation worse in the community.Only when the reinfectivity of the disease is sufficiently weak, the acceptance of visitors may induce the endemicity if the community accepts the visitors with a sufficiently low immune proportion.Furthermore, when the reinfectivity is high, the acceptance of a sufficiently large number of visitors may induce the elimination of the disease if the community can regulate to accept the visitors with a sufficiently high immune proportion.
The visitors certainly play a role of recruitment of hosts for the infectious disease spread in the community.The visitors with a higher susceptible proportion could be regarded as a larger supply of susceptible individuals in the community, and they subsequently provide a fast recruitment of infectives.In contrast, the visitors with a high immune proportion cause only a slow recruitment of infectives with the reinfection.For these reasons, the influence of the visitor acceptance on the epidemic dynamics with a reinfectious disease must depend on the immune proportion in the visitors at the entry.On the other hand, the entry of many visitors could induce a dilution of the infective density in the community at the same time, which is regarded as an advantageous influence of the visitor acceptance against the disease spread.In the epidemic dynamics with our model, a balance of these counteractive As a result, a preferable acceptance of visitors must be regulated to have a sufficiently large immune proportion according to the public health in the community.In this sense, the best policy for the visitor acceptance would be to allow the entry only for the immune visitors.From the results on our model, such an acceptance of only immune visitors may lower the endemic size, and further suppress the endemicity to induce the elimination of the disease spread in the community.
In contrast, when the community was on the way to the disease-eliminated equilibrium before starting the acceptance of visitors, the acceptance of visitors without any epidemiological regulation may cause the revival of the disease spread in the community.Such a case would occur by reduced cautiousness of the disease before starting the acceptance of visitors, which could be caused by the reason that the number of infective residents became rather small in comparison to that at the outbreak.
Our model could be regarded as a consideration on the epidemic dynamics in a season.In this sense, the number of visitors may be beyond the number of residents in the community (i.e.,  > 1 ), as some popular touristic local places like Venice in the vacation season, or a certain place attracting visitors like a newly found gold mine.As another example, we could consider a community accepting many evacuees from a certain calamity.Even though the number of visitors would be smaller than the number of residents in most cases (i.e.,  < 1 ), our results imply that the influence of the visitor acceptance could depend on the infectivity and reinfectivity in the epidemic dynamics, and the regulation on the epidemiological nature of accepted visitors.
Since the infectivity and reinfectivity are not only determined by the nature of disease itself but also by social custom, sanitary condition, and people's behavior (Ferguson et al. 2005;Heymann 2005;Funk et al. 2009;Perisic and Bauch 2009;Manfredi and D'Onofrio 2013;Bavel et al. 2020;Kapitány-Fövény and Sulyok 2020;Thu et al. 2020), the influence of the visitor acceptance could depend also on social factors in the community which accepts the visitors under the epidemic dynamics.Such social factors could be affected by the situation of disease spread during the epidemic dynamics in the community.For example, some strategic/nonstrategic transmission of information about the disease spread or a public health campaign to prevent the further disease spread could alter people's social behavior, and subsequently the risk of infection/reinfection.Hence if the infectivity and reinfectivity would be changeable in the epidemic dynamics, the influence of the visitor acceptance must be qualitatively changed.
The results from our model imply such a possibility that a shift of the infectivity and/ or reinfectivity to the weaker would induce an epidemic situation in which the acceptance of visitors causes the increase in the epidemic size or the revival of disease spread even with the endemicity.If so, there would be repetitive revivals of disease spread in the community, driven by a temporal shift of the infectivity and/or reinfectivity which could bring a feedback influence on the policy to control the disease spread in the community.Such theoretical/mathematical researches on the relation between the disease spread and the nature of hosts are interesting and require further development.

Appendix A Derivation of R 0
In our model, the new cases consist of residents and visitors.Therefore, from the conceptual definition and the mathematical feature of the basic reproduction number, we shall derive it here from the following conditions: for the initial condition (3) with 0 < I v (0) + I r (0) = 0 + I r0 ≪ 1 .This is because the basic reproduction number is defined as the expected number of new infectives produced by one infective individual in an environment consisting of only susceptibles.Following the assumption H3 in Sect.2, we shall adopt the initial condition (3) in order to define the basic reproduction number.Then, in place of the above conditions, we can use the followings: for the initial condition (3) with 0 < I v (0) + I r (0) = 0 + I r0 ≪ 1 .Remark that, in this context about the situation to define the basic reproduction number of our model, the initial infective must be a resident, which matches the assumption H3 in Sect. 2.Then, making use of (3), we have Thus we find that Consequently we can define the basic reproduction number R 0 as follows: This formula can be rewritten as given by (4) to clarify the meaning.
Appendix B Proof for Theorem 6.1 From the arguments in the first paragraph of Sect.6, we find that the dynamics given by ( 5) with = 0 necessarily approaches the dynamics with the following limiting sys- tem in terms of the visitor population: The feasible equilibria are Ẽ0 (1 − , 0) and Ẽ+ (x * v , ỹ * v ) with ( 6).The former Ẽ0 cor- responds to the disease-eliminated equilibrium for the system ( 5), E 00 (1 − , 0, 0, 0) , and so does the latter Ẽ+ to the endemic equilibrium E +0 (x * v , ỹ * v , 0, 0) .The endemic equilibrium E +0 can exist when and only when the condition ( 7) is satisfied.By the local stability analysis with the eigenvalues of the Jacobi matrix at the equilibrium, we can easily find that the endemic equilibrium E +0 is locally asymptotically stable when it exists.In the following part, we shall consider its global stability.
First we set the following mathematical result on the boundedness for the solution of the system (B1): Proof We can obtain the following features from (B1) for the initial condition xv (0), ỹv (0) ∈ D: The first inequality indicates that xv cannot reach 0 from the initial value � x v (0) > 0 .The second equation indicates that ỹv is necessarily positive for any  > 0 and � y v (0) > 0 .Then, from (B1) for � y v > 0 , we can find that xv + ỹv cannot become xv + ỹv = 1 − for any  > 0 and initial condition in D. ◻ Lemma B.1 means that the domain D is invariant for the dynamics given by (B1).Further, when the endemic equilibrium Ẽ+ exists, satisfying the condition (7), we find from ( 6) that Hence we have the following result: exists, it must belong to the domain D defined by (B2).
When the condition (7) is not satisfied, that is, when the endemic equilibrium does not exist, we can easily find that the disease-eliminated equilibrium is globally asymptotically stable, making use of the isocline method shown by Fig. 14(a).In contrast, as seen from Fig. 14(b), when the condition ( 7) is satisfied and the endemic equilibrium Ẽ+ exists, the stability cannot be determined only by the isocline method.
When the endemic equilibrium Ẽ+ exists, let us consider the following function of (x v , ỹv ) in the domain D defined by Lemma B.1: It can be easily found that Further, making use of (B1), we can derive which becomes zero only for xv , ỹv = x * v , ỹ * v , and negative for any xv , ỹv ≠ x * v , ỹ * v in D. These features of V indicates that it is a Lyapunov function according to the endemic equilibrium Ẽ+ for the system (B1).Therefore, when the endemic equilibrium Ẽ+ exists for the system (B1), it is globally asymptotically sta- ble with respect to the dynamics given by (B1) with the initial condition in D. Since the dynamics of ( 5) with = 0 must eventually approach that of (B1), this result shows the global stability of the endemic equilibrium for (5).

Appendix C Proof for Theorem 8.1
The characteristic equation det(J 00 − E) = 0 with the Jacobi matrix J 00 accord- ing to the disease-eliminated equilibrium E 00 (1 − , 0, 0, 0) for the system (5) becomes Hence we have two negative eigenvalues −1 and − with the other two given by the roots of the equation that the above 2 × 2 determinant is equal to zero.Both of them have negative real parts if and only if that is, Therefore we find that the second inequality gives the necessary and sufficient condition that every eigenvalue has a negative real part.It can be expressed as the inverse inequality of (9).At the same time, we can find that, if the inverse of the above second inequality is satisfied, there exists an eigenvalue with a positive real part.Then the disease-eliminated equilibrium E 00 is unstable.

Appendix F Proof for 8.3
The Jacobi matrix about the endemic equilibrium E ++ becomes where Then the characteristic polynomial for the eigenvalue about E ++ can be obtained as  given by (F6) to determine the other three eigenvalues for E ++ .Every coefficient of h( ) is positive: a 2 > 0 , a 1 > 0 , and a 0 > 0 .Since 0 < Φ < 1 , we have and subsequently find that Consequently from the Routh-Hurwitz criterion, we can find that all roots of h( ) = 0 have negative real parts.Therefore it has been proved that every eigenvalue for E ++ has negative real part.This result shows Theorem 8.3 about the local stabil- ity of E ++ when it exists.

Now we can derive
Thus we can result that ( c ) < 0 if and only if  <  c where c is defined by ( 19).That is,  * ∕ > 0 if and only if  <  c .Lastly we have Lemma G.3 x * v > 0 is monotonically decreasing in terms of  > 0 if and only if  <  c when  < 1 and R 00 > 1.
It can be easily found that  c < 1 when R 00 > 1.
On the other hand, from (10), we can easily find that the endemic sizes y * r , y * v , and z * are monotonically decreasing in terms of x * v > 0 .Therefore, from Lemma G.2, we can get the following lemma: Lemma G.4 The endemic sizes y * r , y * v , and z * are monotonic in terms of when  < 1.

Fig. 1
Fig.1Scheme of the model for the epidemic dynamics in a community accepting temporal visitors, given by the system of (1) and (2)

Fig. 10 3
Fig. 10 Parameter region and boundary indicated by the condition the condition (9) with the results given by Theorems 8.1, 8.2, and 8.3 when all visitors accepted by the community is susceptible with = 0 : a R 00 > 1 + c ; b R 00 = 1 + c ; c R 00 < 1 + c .Numerically drawn with a c = 1.0 ; b c = 3.0 ; c c = 5.0 , and commonly R 00 = 4.0 .Refer to Sect.9.3

Fig. 13 3
Fig. 13Classification of the parameter region of ( , ) according to the -dependence of the change in the endemic size.Numerically drawn with R 00 = 4.0 and c = 1.0 .Regions Ω ± correspond to those in Fig.6

Fig. 14
Fig.14Application of the isocline method for the system (B1) when the condition (7) is (a) not satisfied; (b) satisfied